Soliton solution of the osmosis K(2, 2) equation 
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Abstract 

In this Letter, by using the bifurcation method of planar dynamical systems, we 

in ' 

PLh ' obtain the analytic expressions of soliton solution of the osmosis K(2, 2) equation: 



u t + (u 2 ) x - (u 2 ) xxx = 0. 
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Since the theory of solitons has very wide applications in fluid dynam- 
ics, nonlinear optics, biochemistry, microbiology, geophysics and many other 
fields, the study of soliton solutions has become one of the important issues 
of nonlinear evolution equations [l]-[8]. 

Recently, Xu and Tian [9] investigated the osmosis K(2, 2) equation 

u t + (u 2 ) x - (u 2 ) xxx = 0, (1.1) 
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where the positive convection term (u 2 ) x means the convection moves along 
the motion direction, and the negative dispersive term (u 2 ) xxx denotes the 
contracting dispersion. They obtained the peaked solitary wave solution and 
the periodic cusp wave solution of Eq.(l.l). Unfortunately, the results are not 
complete. In the present Letter, we shall continue their work and obtain the 
smooth soliton solutions of Eq.(l.l), so that we can supplement the results of 



The remainder of this Letter is organized as follows. In Section 2, we state 
the main results which are analytic expressions of the smooth soliton solutions 
of Eq.(l.l). In Section 3, we give the proof of the main results. 

2 Main results 

We state our main result as follows. 
Theorem 2.1 For given constants c ^ and — y < g < — ^j-, let £ = x — ct, 
(1) If c < 0, then Eq.(l.l) has a soliton solution of the form: 



[9]. 
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where 



fh(<p) 
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We shall give the proof of this theorem in Section 3. Now we take a set of 
data and employ Maple to display the graphs of the soliton solutions oem-( 
plane. 



Example 2.1 Taking c = —2 and g = —0.9999 (corresponding to (1) of 
Theorem 2.1), it follows that ip^ = -1.01, ip x = -0.979924, h = 0.646667 
and l 2 = -0.326567, l 3 = -1.30627, a x = 0.0404, b x = -1.37333, and a x = 
—5.02494. We present the graph of the soliton solution in Fig.l (a). 

Example 2.2 Taking c = 2 and g = —0.9999 (corresponding to (2) of The- 
orem 2.1), it follows that <p$ = 1.01, (pf = 0.979924, m x = -0.646667 and 
m 2 = -0.326567, m 3 = -1.30627, a 2 = 0.0404, b 2 = 1.37333, a 2 = 5.02494. 
The graph of the soliton solution is presented in Fig. 1(b). 




(a) (b) 

Fig. 1. The graphs of the soliton solutions, (a) c = —2, g = —0.9999; (b) c = 2, 
g = -0.9999. 



3 Proof of main results 

Eq.(l.l) also takes the form 

u t + 2uu x — Qu x u xx — 2uu xxx = 0, (3-1) 
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Let u = </?(£) with £ = x — ct(c ^ 0) be the solution of Eq.(3.1), then it 
follows that 



- ap' + 2iptp' - Gipip" - 2tp<p'" = 0. (3.2) 

Integrating (3.2) once we have 

-ap+((p) 2 -2(<p') 2 -2<p<p" = 9, (3.3) 

where g is the integral constant. 

Let y = ip', then we get the following planar dynamical system: 



^ = v 



dy _ (p 2 -ap-g-2y 2 



(3.4) 



with a first integral 

H(<p, y) = V 2 {y 2 - \p 2 + ~<p+ \g) = h, (3.5) 



where h is a constant. 

Note that (3.4) has a singular line ip = 0, to avoid the line temporarily we 
make transformation d£ = 2<pdQ. Under this transformation, Eq.(3.4) becomes 



% = *pv 

(3.6) 

^ = cp 2_ C(p _ g _ 2y 2 



Eq.(3.4) and Eq.(3.6) have the same first integral as (3.5). Consequently, sys- 
tem (3.4) has the same topological phase portraits as system (3.6) except for 
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the straight line p = 0. Obviously, ip = is an invariant straight-line solution 
of system (3.6). 

For a fixed h, (3.5) determines a set of invariant curves of (3.6). As h is 
varied, (3.5) determines different families of orbits of (3.6) having different 
dynamical behaviors. Let M(ip e , y e ) be the coefficient matrix of the linearized 
system of (3.6) at the equilibrium point (<p e ,y e ), then 



\ 



Ve 2y3 e 
2<p e - c -Ay e 



and at this equilibrium point, we have 



J(<p e ,y e ) = detM(<p e ,y e ) = -Ay 2 e - Ap> e (ip e - -) 



p(<Pe,ye) = trace(M(^ e ,y e )) = -3y e 



By the theory of planar dynamical system (see [10]), for an equilibrium point 
of a planar dynamical system, if J < 0, then this equilibrium point is a saddle 
point; it is a center point if J > and p = 0; if J = and the Poincare index 
of the equilibrium point is 0, then it is a cusp. 

In [9], Xu and Tian reported that when — ^ < g < — system (3.4) has 
two equilibrium points ((po, 0) and (ipg , 0). (i) If c < 0, then (<fo, 0) is a saddle 
point while (<Po,0) is a center point. There is inequality <Pq < | < <p>Q < 0; 
(ii) If c > , then (</?o,0) is a center point while (v?o,0) is a saddle point. 
There is inequality < <p^ < | < ip$ . 

In the parameter region: — y < g < — ^j-, we show the phase portraits of 
system (3.4) in Fig. 2. 



6 



(a) (b) 

Fig. 2. The hase portraits of system (3.4) when — ^ < g < — ^j-. (a) c < 0; (b) 
c> 0. 



Usually, a soliton solution of Eq.(l.l) corresponds to a homoclinic orbit of 
system (3.4). In Fig. 2 (a), the homoclinic orbit of system (3.4) can be expressed 
as 



2/ = ± ^ /or <A) < V < <Pi , (3-7) 



where </?q , ip±, l\ and ^ are in (2.3), (2.4), (2.5) and (2.6), respectively. Sub- 
stituting Eq.(3.7) into the first equation of system (3.4) and integrating along 
the homoclinic orbits, we have the expression of the solion solution as in (2.1). 

In Fig. 2(b), the homoclinic orbit of system (3.4) can be expressed as 



V = ± ^ for if{ < <p < <p%, (3. 



where <Pq, <pf, ui\ and m 2 are in (2.13), (2.14), (2.15) and (2.16), respectively. 
Substituting Eq.(3.7) into the first equation of system (3.4) and integrating 
along the homoclinic orbits, we have the expression of the soliton solution as 
in (2.11). The proof of Theorem 2.1 is completed. 
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